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1. Introduction 

Feynman p.] famously foreshadowed atomic-scale machines that could perform similarly 
to their macroscopic analogs. Nowadays the problem of designing nanomechanical 
devices, in particular, to reduce friction by means of nano- and micro-bearings [2] is 
real. The possible use of fullerene Ceo for molecular nanobearings gave rise to Molecular 
Dynamics (MD) modeling IH [3J E]. Other fullerene-like metal dichalcogenide MX 2 
(where M=Mo or W and X=S or Se) molecules were considered as additives in oil 
lubricants, and predicted to provide interesting tribological properties (e.g., see [7] and 
references therein). Simulations showed that ball-shaped molecules may either slide or 
rotate over a surface, depending on the substrate and the position of the molecule. For 
example, Cgo slides on graphite in the AB configuration (hexagonal Cqq ring lying flat on 
graphite), but rotates in the on-top frustrated AB configuration where one Cgo corner 
atom faces the center of the graphite hexagon. The rolling configuration is characterized 
by very low friction 013], with a predicted friction coefficient of the order fi ~ 0.01 — 0.02 
[I] or even smaller 0. 
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Attempts to realize these ideas experimentally had only limited success so far. A 
single Cgo molecule confined between two solid substrates may begin to roll when a 
torque of order 1CT 19 Nm is applied [8]. However, Ceo molecules actually condense to 
form close-packed layers, as found, e.g., on a graphite substrate 0, [TP] . A single C 60 
monolayer (ML) takes a crystalline structure with 2D spatial order at low temperatures. 
It undergoes a first-order orientational order- disorder transition at T = T m f« 260 K [llj, 
the molecules exhibiting free rotation at T > T m . At T m there is an abrupt change in 
friction [12], but the lowest friction coefficient is of order \i ~ 0.15 [TQl H2], worse than 
with traditional oil-based lubricants. Coffey and Krim [13J reported a quartz crystal 
microbalance study of one or two Ceo monolayers adsorbed on Ag(lll) or Cu(lll). 
There are no rotations in a Cgo ML on Cu(lll), and only a slow change of molecular 
orientations in the C6o/Ag(lll) ML. For two ML instead, Ceo molecules in the second 
layer rotate freely at 300 K. However a molecularly thin methanol film deposited over the 
C 60 , failed to show either the expected low friction, or any essential difference between 
these systems. Thus this particular nanobearing design apparently would not work. 

Some charge transfer and bonding between Cgo and the metal substrate may be held 
responsible for hindering the rolling. Another reason lies in the full layer coverage. Balls 
in macroscopic bearings are arranged so as to prevent contact, but rolling molecules in 
the ML are always in contact, hindering their mutual rolling and jamming the same 
way two ingrained rolling cogwheels would. As discussed earlier on [Hj , a way to avoid 
jamming is to lower dramatically the coverage, well below one ML (the molecule density 
should anyway be sufficient to prevent the two surfaces from touching). In view of that, 
and in the lack of well defined low coverage experiments, a study of the single molecule 
rolling friction represents a natural starting point, and indeed a revealing one. 

Macroscopically, the main source of rolling friction of a ball or tire comes from 
deformation. Both substrate and roller are (elastically or plastically) deformed at the 
contact. The deformation energy is partly released and lost as bulk frictional heat when 
the roller moves on |15j . By designing the bulk so that dissipation is poor, rolling 
friction can be made 10 2 to 10 3 times lower than the sliding friction; the latter being 
due to adhesion, i.e., breaking and re-forming of slider-substrate bonds. It the previous 
work [14] we studied molecular rolling friction for the system, where the lubricant and 
both substrates were constructed of the same molecules, so that the lubricant and 
substrates were deformable and commensurate. The minimal friction coefficient found 
in simulation, was of order fi > 0.15 in agreement with available experimental data. 
Naturally it emerges a question, what would be a value of /i for the rigid substrates, 
when the losses due to deformation are absent, or for the case of incommensurate 
lubricant /substrate interface. As the roller size is decreased however, adhesion grows in 
importance, eventually becoming the main source of friction. To rotate a molecule, one 
has to break the molecule-substrate bonds from one side of the molecule and create new 
bonds on the opposite side. Thus, there are no reasons to expect that molecular rolling 
friction should be much lower than the sliding friction. 

Our present goal is to understand what could be the lowest friction coefficient 
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attainable for molecular rolling and which system parameters might provide it. Besides, 
we show that a macro- to microworld mapping does work, but one has to choose properly 
the macroscopic counterpart, which in the present case is a cogwheel. Because we are 
interested in general trends, we explore a minimal two-dimensional (2D) model, which 
allows us to span a large number of parameters, and also provides an easier visualization 
of the processes inside the lubricant. 



2. Model 



We consider two substrates with lubricant molecules in between, all of them made up 
of classical point particles (atoms). Atoms can move in the (x,y) plane, where x is 
the sliding direction and y is perpendicular to the substrates. The substrates, pressed 
together by a load force Fi = N s fi, consist of rigid atomic chains of length N s and 
equal lattice constant R s , so that the system size in the sliding direction is L x = N S R S 
and the total mass of the substrate is N s m s (we use periodic boundary condition along 
x). The bottom rigid substrate is fixed at x = y = 0, the top one is free to move in 
both x and y. The top substrate is driven along x with speed v s through a spring of 
elastic constant k s . The spring force F, whose maximum value before motion measures 
the static friction force F s , and whose average during smooth motion F^ = (F) is 
the kinetic friction force, is monitored during simulation (throughout the paper we 
normalize forces per substrate atom / = F/N s ). Thus, our model is a 2D variant of a 
typical experimental setup in tribology [15], [16]. Between the substrates we have circular 
("spherical") lubricant molecules built as in Ref. [Hj. Each molecule has one central 
atom and L atoms on circle of radius R m = Rn/2sm(ir/L) so that their chord distance 
is Ru. They are coupled with the central atom, additionally to the 12-6 Lennard- Jones 
(LJ) potential, by stiff springs of elastic constant K m , V sta h(r) = \K m {r — Rstab) 2 , where 
the distance R stah = R m + (12 V u /K m R m ) [(R n /R m f - (R a /R m ) 12 ] is chosen so that the 
total potential Vlj(^) + V r sta b(^) is minimum at r = R m . With K m = 100 the resulting 
stiff molecular shape resisted destruction during the simulations. All atoms interact via 
the LJ potential Vlj(?~) = V aa > {R a a'/ r ) 12 — 2 {R a a'/ r ) 6 , where a, a' = s or I for the 
substrate or lubricant atoms respectively. Thus, the lubricant-lubricant interaction is 
described by the parameters Vu and Ru, while the lubricant-substrate interaction, by V s i 
and R s i (direct interaction between the top and bottom substrates is omitted, as they 
are not allowed to touch). We use dimensionless units, where m s = mi = 1, Ru = 1, 
and the energy parameters V aa > takes values around V aa i ~ 1. 

Because a 2D model cannot reproduce even qualitatively the phonon spectrum of 
a 3D system, and because frictional kinetics is generally diffusional rather than inertial, 
we use Langevin equations of motion with Gaussian random forces corresponding to 
temperature T, and a damping force f VjX = —m r)(y) x — m r](Y — y){x — X), where x, y 
are the atomic coordinates and X, Y are the coordinates of the top substrate (the force 
f ViV is defined in the same way) . The viscous damping coefficient is assumed to decrease 
with the distance from the corresponding substrate, rj(y) = 770 [1 — tanh(y/yd)], where 
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typically 770 = 1 and y d ~ 1. 

We present simulation results for molecule friction from L = 5 (the simplest circular 
molecule) up to L = 13 and 14, which may be considered as a 2D version of fullerenes. 
In fact in the 3D case, the surface area of the spherical molecule is s = AnR^. If we 
put L 3 = 60 atoms on the surface, this gives s ~ L^R U) or R m /Ru « 2.18. In 2D, the 
length of the circle is 2nR m « LR U) or L « 2nR m /Ru, which leads to L « 13.7 for the 
same ratio R m /Ru as for 3D fullerenes. 



3. Rigid molecule 

We first consider a rigid circular molecule, i.e., Vu = 00 and K m = 00. Let us fix X of 
the top substrate and seek minimum of the potential energy V by varying the coordinate 
Y of the top substrate, and the center (x c , y c ) and the rotation angle of the molecule. 
The X dependence of V, (x c , y c ), and defines the adiabatic trajectory, which describes 
the joint substrate and lubricant motion when infinitely slow. We define the activation 
energy E a = max [V(X)] — min [V(X)], and the magnitude of the static friction force, 
approximated as f s = max [dV(X)/dX] (f s ~ E a in our units). 
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Figure 1. Activation energy E a as a function of the ratio of the substrate lattice 
constant R s to Ra for the rigid L = 6 molecule, for /; = 0.5, V s i = 1/9, and R s i = Ru. 
Open symbols correspond to sliding, solid symbols to rolling. 



Figures [T] and [2] show the results for the L = 6 molecule when, to simplify further, 
R s i is kept constant, R s i = Ru. The energy V(X) is periodic with R s (or a multiplier 
of R s ). The molecular angle 4> varies by A0 as the potential energy V(X) changes from 
minimum to maximum. Because 4>(X) must be continuous, the motion corresponds to 
sliding if A<p < </)q = 2n/L, while if A<ft > <po the molecule must rotate when it moves. 
As figure [1] shows, for R s < Ru the motion corresponds to sliding, i.e., the molecule is 
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Figure 2. Sliding adiabatic motion of the rigid L = 6 molecule for R s /Ru = 0.66 
(A<f> < 4>o, left panel) and rolling for R s /Ru = 1.29 (A</> > </> , right panel). Other 
parameters as in figure [TJ Lower panels: X-dependence of potential energy V(X) 
and the rotation angle <p(X)/<po, where 4>q = 2it/L. Top panel: configurations as the 
molecule moves from one minimum of V(X) to the next. 

shifted as a whole, slightly oscillating during motion (figure [2], left panel). Similarly to 
the motion of a dimer in a periodic potential [32], the activation energy has maxima at 
Ru = nR s (where n is an integer) and minima at Ru — (n— 1/2)R S . On the other hand, 
for Ru < R s the motion corresponds to rolling (figure El right panel). Here E a (R s ) has 
minima at some values of the ratio R s /Ru (e.g., for R s /Ru ~ 1.29 in figured]). 

Varying R s in figure (TJ we kept fixed the equilibrium distance R s i for the lubricant- 
substrate interaction. More realistically, it might be reasonable to set, e.g., R s i = R s , 
in which case, as we observed, the interval of R s values where rolling prevails is wider 
than for fixed R s i. Further preference for rolling over sliding is found for increasing load 
fi and for decreasing interaction strength V s i. We also note that when sliding wins over 
rolling for R s < Ru, it provides a lower activation energy. Recalling that (fio = 2n/L, the 
region of parameters for rolling should increase with L - a rounder wheel rolls better. 
The R s dependence of E a (R s ) for increasing size L (figure [3]) shows rolling for all R s and 
for all L > 5, except for L = 6 which shows both rolling and sliding (see open symbols in 
figure [3^). As R s varies, the value of E a changes by more than two orders of magnitude 
for even L and more than three for odd L, with deep sharp minima separated by broad 
maxima. Clearly, by suitably choosing R s /Ru a very strong decrease of rolling friction 
is attainable. 

The deep minima of E a (R s ) are explained by simple engineering - a "cogwheel 
model". Consider the molecule as a cogwheel with L cogs, primitive radius R m and 
external radius R* = R m + h, where h oc R s i. The chord distance between nearest 
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Figure 3. Rigid molecule activation energy E a versus R s /Ru for /; = 0.5 and 
V s i — 1/9. Unlike figure [IJ here R s i — R s . (a) is for even L = 6, 8, 10, 12 and 
14; (b) is for odd L = 5, 7, 9, 11 and 13. Empty triangles in L=6 indicate sliding 
motion intervals. 

cogs is = 2R*sm(n/L). Best rolling conditions are expected when R* t matches 
the substrate potential period R s , i.e., for R^ = R^ and its fractions, R^ = R*J2, 
i?(3) = i?*/3 3 etc. The main minimum of E a (R s ) is expected at 

R^/Ru = 1 + (2h/R u ) sin(7r/L). (1) 

As shown in figure HI the cogwheel model (OQ) with h = /3R s i, where (3 is a parameter, 
can fit very well the shift of minimum position with molecular size L. It can explain 
its variation with load (the radius R* and therefore h decrease as the load grows) as 
well as with the lubricant-substrate interaction V s i (R* and h decrease with V s i). It also 
accounts for the even-odd effect since odd L involves ingraining perfectly one substrate 
at a time, justifying why roughly double values of j3 are needed for even relative to 
odd L. 
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Figure 4. Position Ri of the main minimum of E a (R s ), extracted from figure [3l for 
increasing molecular size L. Curves are fits to the cogwheel model (flj. The asymptotic 
limit of 1 is still relatively far for reasonable molecular radii. 

4. MD simulation 

The simulation results for the static friction of a deformable circular molecule are 
presented in figure El As one could expect, the L = 3 or 4 "circular" molecule does 
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Figure 5. The static friction coefficient /i s — f s /fi and the friction force f s (inset) 
for a single circular molecule as a function of size L. The parameters are V s i — 1/9, 
R si = R s , fi = 0.5, N s = 19, R s = 2/3, v s = 0.03, k s = 10~ 3 , 770 = 1, and y d = R s . 

not roll; instead we observed its "creep" with a relatively large friction. For larger 
values of L, L > 5, the molecule may either roll or slide. For rolling in the case of even 
values of L (L — 6, 8, 10, 12 and 14) one needs to break simultaneously two lubricant- 
substrate bonds (one connecting the lubricant molecule with the bottom substrate, and 
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Figure 6. The static friction force f s for the deformable circular molecule (left axes, 
K m = 100 and Vu = 1, solid curve and circles for v s = 0.003 and stars for v s = 0.0003) 
and the activation energy E a for the rigid molecule (right axes, K m — oo and Vu = oo, 
red open symbols and dotted curve) as functions of R s /Ru for /; = 0.5, V s i — 1/9 and 
R s i = R s . (a) is for L = 14, and (b), for L = 13. 



one, with the top substrate). Therefore, f s should be approximately independent of L, 
as indeed is observed in simulation for L > 8. For odd values of L, L = 5, 7, 9, 11, 13 
and 15, f s is at least two times smaller than for a nearest even L value, because one 
needs to break one bond only at a time. For all L > 7 the static friction is relatively 
low, /j, s < 0.1, and for large odd values the friction may reach quite low values. 

The results obtained for the rigid molecules in section [3j are qualitatively confirmed 
by the static friction force obtained from simulation with deformable molecules. Figure [6] 
compares the results obtained for the rigid molecule with the MD calculation of the static 
friction force of the deformable molecule. The agreement between these two dependences 
is reasonable, at least qualitatively. The friction coefficient /i s = f s /fi ranges from 
H s ~ 0.1 at R s /Ru ~ 0.7 to fi s ~ 0.01 or even /i s ~ 0.001 at R s /Ru ~ 1.1. These 
results are robust to a change of model parameters. For example, figure [7] compares the 
dependences f s (R s ) for two values of the amplitude of lubricant-substrate interaction, 
V s i = 1/9 and 1/3, and for two values of the load, = 0.5 and 0.1. The next two 
figures show the dependence of the static and kinetic friction on V s i (figure ED and on 
the load (figure [9]) ; the latter demonstrates that the friction force approximately follows 
the Amontons law 

fs,k ~ fosflk + ^s,kfl- (2) 

Visualization of MD trajectories shows that for R s /Ru = 0.7, where friction is high, 
rolling rotation is accompanied by a molecular shift/sliding, - much as cogwheels with 
excessive clearance would do - while for R s /Ru = 1.1, where friction is low, the motion 
corresponds to pure rotation, corresponding to optimal cogwheel coupling. 

Simulations showed that the results presented above remain valid at nonzero 
temperature T. When T increases, we observed both static and kinetic friction force 
to decrease, the stick-slip changing to creep and finally to smooth motion at a high 
temperature. Moreover, we found a transition from stick-slip to smooth rolling for 
increasing velocity (figure fTUl) . The cogwheel effect remains, and for example calculated 
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Figure 7. The static friction force f s as a function of the substrate lattice constant 
R s for (a) L = 14 and (b) L = 13 for different system parameters: (i) /; = 0.5 and 
Vsi = 1/9 (solid curve and circles and stars), (ii) fi = 0.5 and V s i = 1/3 (down triangles 
and red dotted curve and stars), and (Hi) ft = 0.1 and V s i — 1/9 (up triangles and 
blue dotted curve). Other parameters are R s i = R s , K m = 100, and Vu = 1. 



static friction for R s /Ru = 0.7 and R s /Ru = 1.1 still differ by a factor of 10 or more. 
The critical velocity v c of the transition from stick-slip to smooth rolling also differs by 
a factor of about four in the two cases. Moreover we always find fk ^ fs- 

The present approach to the single rolling molecule can be extended to a finite 
coverage of lubricant molecules. For example, figure [TT] shows the friction force for a 
finite concentration of lubricant molecules, which may be compared with those of figure [5] 
for a single molecule. These results are for approximately the same load per one lubricant 
molecule (fiN s /M w 9.5 in both cases), and we used a relatively low concentration of 
lubricant molecules, M/N s rs 0.05, to avoid jams. The dependence f s (L) in figure [Til 
is essentially similar to that of figure El although the even-odd oscillation of f s with 
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Figure 8. The static friction force f s (solid curves and symbols, stick-slip motion) 
and the kinetic friction force fk (dotted curves and open symbols, smooth rolling) as 
functions of the amplitude V s i of lubricant-substrate interaction for the deformable 
L = 14 circular molecule for two values of the ratio R s /Ru = 0.7 (down triangles) and 
Rs/Rn = 1-15 (blue up triangles). The parameters are fi = 0.5, N s = 19, k s — 10 -3 , 
Vo = 1, Vd = R s , Rsi = R s , K m = 100, and V u = 1. 
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Figure 9. The static friction force f s (solid curves and symbols, stick-slip motion) 
and the kinetic friction force fk (dotted curves and open symbols, smooth rolling) as 
functions of the load /; for the deformable L = 14 circular molecule for two values of 
the ratio R s /Ru — 0.7 (down triangles) and R s /Ru — 1.15 (blue up triangles and red 
stars). The parameters are V s i = 1/9, N s = 19, k s = 10~ 3 , 770 = 1, yd, = Rs, Rsi = Rs, 
K m = 100, and V u = 1. 



L are less pronounced at the finite concentration because of collisions between the 
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Figure 10. The friction force as a function of the driving velocity v s for two values 
of the ratio R s /Ru = 0.7 (down triangles) and R s /Ru = 1.1 (blue up triangles). Solid 
symbols are for the static friction, open symbols correspond to the kinetic friction; v c is 
the critical velocity of the transition from stick-slip to smooth rolling. The parameters 
are L = 14, /, = 0.5, V sl = 1/9, N s = 19, k s = lO" 3 , Vo = 1, y d = R„ R,i = R,, 
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Figure 11. The static (diamonds) and kinetic (red circles) friction force for a finite 
concentration of circular molecules. The parameters are N s — 151, /; = 0.5, M = 8, 
v s = 0.03 for the static friction (stick-slip) and v s — 0.3 for the kinetic friction (smooth 
motion); other parameters are as in figure Small blue circles and dotted curve show 
the results for the single molecule from figure [5j 



molecules. As for kinetic friction at high driving velocity v s = 0.3 for smooth motion, 
it demonstrates a more monotonic behavior with L without even-odd oscillations. The 
function fk{L) reaches a minimum at L = 6 where < 0.01, and then increases until 
L = 12; at higher values of L the dependence fk{L) approximately repeats the behavior 
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Figure 12. Dependence of the kinetic friction force at v s — 0.3 on the number 
of circular lubricant molecules M (bottom axes) or on the dimensionless coverage 
Q M = M/Mi (top axes). The parameters are: L = 14, N s = 151 ) = Q.2, 



R s /R u = 1.15, V sl 
and T = 0.05. 



1/9, R sl = R s , k s 



io- 



1, y d = R s , K n 



100, 



Finally, figure [12] shows the dependence of the friction force on the concentration 
of lubricant molecules for L = 14 and R s /Ru = 1.15, which provided a low friction 
in the single molecule case (figure [7^). When M increases, the total loading force 
Fi = fiN s is split over the M molecules, so that for a given molecule the load is 
fn — fi/M. As the load decreases with M, the friction force per molecule f s i,ki should 
also decrease according to (j2J). At the same time, the total friction force should increase, 
f S) k = Mf s i t ki- A combined effect is a slow increase of the friction with M as shown 
in figure [12] with dotted curve and open symbols. In a real situation, coalescence may 
lead to jamming, with molecules blocking their mutual rotation [14]. In our model, 
jamming starts already at 9m ~ 0.1 and completely destroys rolling at 9m > 0.3 (here 
9m = M/Mi is the coverage, with Mi the number of molecules in the monolayer). 



5. Conclusion 



Summarizing, we can extract from our 2D model the following conclusions. Rolling 
spherical lubricant molecules can indeed provide better tribological parameters than 
sliding atomic lubricants. The effect may be as large as in macroscopic friction, where 
rolling reduces friction by a factor of 10 2 — 10 3 , however only for sufficiently low 
coverage of lubricant molecules, and for specially chosen values of the ratio R s /Ru, 
corresponding to perfect cogwheel rolling. To check experimentally these predictions, 
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it would be interesting to study friction coefficient for different spherical molecules, 
different coverages, and different substrates. Also, the relative ingraining between the 
rolling molecule and the substrate may be improved by adjusting the applied load, 
as it was demonstrated experimentally for the molecular rack-and-pinion device [18J. 
Inert nonmetal surface (such as perhaps self-assembled monolayers) may represent a 
better choice of substrate than metals for fullerenes deposition. Because 3D rolling 
has an azimuthal degree of freedom, the cogwheel effect described should be direction 
dependent, and rolling friction should exhibit anisotropy depending on direction. 
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